Quantum Limits of Stochastic Cooling of a Bosonic Gas 
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The quantum limits of stochastic cooling of trapped atoms are studied. The energy subtraction due to the 
applied feedback is shown to contain an additional noise term due to atom-number fluctuations in the feedback 
region. This novel effect is shown to dominate the cooling efficiency near the condensation point. Furthermore, 
we show first results that indicate that Bose-Einstein condensation could be reached via stochastic cooling. 
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Ultracold atomic gases are generated by laser cooling (J], 
evaporative cooling |2j] or sympathetic cooling |3]. These 
techniques have been successful in preparing Bose-Einstein 
condensed |4] and Fermi degenerate |5] atomic gases. These 
states of matter have proven important in both the fundamen- 
tal physics of weakly interacting many-body systems 0] and 
in applications |7|,|S|]. Nonetheless, these cooling methods 
have some limitations. For example, laser cooling requires 
closed-cycle transitions, evaporative cooling leads to a loss of 
a significant fraction of the atoms, and sympathetic cooling 
requires careful selection of buffer species with sufficiently 
large scattering cross sections. 

A much more general strategy that avoids these limitations 
is stochastic cooling |9]. It is a Maxwell-demon strategy that 
uses information obtained from measurement to coherently 
reduce the energy of part of the system. Stochastic cooling 
is based on the repeated application of a feedback loop. The 
cooling is obtained by the combination of measurement and 
controlled Hamiltonian interaction during the feedback oper- 
ation. The interaction provides an energy exchange with an 
external field and the preceeding measurement ensures the ir- 
reversibility of this exchange. In this way the feedback mech- 
anism may be thought of as acting as a dissipative reservoir. 
Classically cooling occurs due to the extraction of information 
on the phase-space localization of particles, and its subsequent 
use to reposition the particles, leading to a phase-space com- 
pression 

In high-energy physics it has been employed for cooling 
the transverse degree of freedom of a particle beam A 
measurement of the transverse momentum of a fraction of the 
particles is made. A control field sets the momentum to zero, 
which together with a subsequent remixing of the particles 
leads to phase-space compression and cooling of the trans- 
verse motion. Recently, stochastic cooling was proposed for 
trapped atoms and it was shown that both momentum mea- 
surement and shift could be realized by optical fields 
The required remixing of atoms is provided here by the os- 
cillation in the trap. Moreover, interactions between atoms, 
such as collisions, may provide a further enhancement of this 
remixing. Classical calculations for a ID atomic gas showed 
a pronounced cooling effect I Kill , so that stochastic cooling 
may perhaps be an alternative to standard cooling methods. 



However, to best of our knowledge, it is not yet known 
to what temperatures such a method eventually will cool the 
atoms and whether, for example, Bose-Einstein condensa- 
tion can be reached. Technical heating effects, that are in- 
herent in the proposed optical implementation, have been dis- 
cussed flOfl . However, the fundamental limits of such a cool- 
ing method, due to the discreteness of atoms and their quan- 
tum correlations, seem to be unexplored. A reason for the 
lack of knowledge of these limits may be that many-atom cor- 
relations play a central role. The atoms cannot be treated as 
individual entities, which would considerably simplify a theo- 
retical description, but require a quantum many-atom descrip- 
tion full . 

Measurement in quantum mechanics always leads to a 
back-action in the conjugate variable, and one might expect 
that this will saturate the cooling at ultralow temperatures. It 
will be shown in this Letter that a further fundamental heat- 
ing mechanism arises from the quantum fluctuations of the 
number of atoms in the feedback region. To reach ultralow 
temperatures, this heating has to be circumvented, which will 
be shown to be possible by the choice of the feedback region. 
We treat stochastic cooling in a fully quantum-field theoreti- 
cal way and show initial results that indicate the possibility of 
reaching Bose-Einstein condensation. 

The feedback loop of stochastic cooling consists of a mea- 
surement of the total momentum P w of the atoms in the feed- 
back region and a subsequent interaction in the same spatial 
region that compensates for the observed value P, i.e., shifts 
the total momentum to zero. The spatial region where atoms 
are subject to the feedback is defined by the beam of a laser 
field that implements the measurement and subsequent shift, 
cf. Ref. fioll . For simplicity we consider the case where the 
laser beam is aligned along the z axis and the beam waist w 
shall be a step-like function in the xy-plane, cf. Fig.^ Note 
that in z-direction no spatial restriction is assumed. 

The measurement and subsequent shift shall occur without 
time delay and on timescales much smaller than the intrinsic 
timescale of the system 1 / ujo, where ljq is the trap frequency 
of the 3D isotropic, harmonic potential. Then the entire feed- 
back loop approximately acts instantaneously in time. Fur- 
thermore, we take into account the resolution a of the mea- 
surement of momentum P w , that is determined by external 
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FIG. 1: Geometry of the feedback setup. The laser beam is aligned 
along the z axis, it determines the size of the feedback region. 



optical fields implementing the measurement. 

The feedback acts on the many-atom density operator g as 

g+ = jdPU{P)M{P-P w )Q-M\P-P w )U\P), (1) 

where g T are the many-atom density operators before (— ) 
and after (+) the feedback. The measurement is described 
by the positive operator-valued measure M, where |M(P)| 2 
is a Gaussian of width a. The unitary operator U(P) — 
exp[iPQ w ] implements the subsequent shift of the measured 
momentum P to zero lll2ll . The final density operator is found 
by averaging over all possible measurement outcomes. 

Using the bosonic atom-field operator <fi(r) with commuta- 
tor [$(r), </>^(r')] = <5( 3 )(r— r'), the measured total momentum 
P w and the center-of-mass Q w of atoms in the feedback re- 
gion are 



P w = dA dztf{r)(-id z )<t>(r), (2) 

J W J 

Q w = JdA Jdzft(r)z$(r), (3) 

where dA is integrated over the beam waist w. Their commu- 
tator is [Q w , P w ]=iN w /N e , where the operator of the number 
of atoms in the feedback region is 



N w = dA dz^(r)$(r). 



(4) 



Note that in Eq. Q the estimated number of atoms N c is used 
and not the proper operator N w : Since, roughly speaking, in 
the Schrodinger picture U(P) = Yli ex P{iP/N e qi), where 
qi is the coordinate of the ith atom, each atom is shifted by 
~P/N C . That is, since the true atom number is unknown — it 
is not measured — an estimate for the atom number, N e , has 
to be employed for properly shifting each atom's momentum. 

For cooling we are interested in the average energy that 
is subtracted by the feedback from the set of non-interacting 
atoms, 



AE = (H) 



(H) 



(5) 



Here H is the Hamiltonian of the atoms in the harmonic trap 
potential and (. . )± = Tr(p+ . . .) . After a detailed calcula- 
tion, using the results of Ref. | ljj], we find that AE consists 
of both cooling and heating terms. 

Firstly, consider an expansion of the kinetic energy of the 
center of mass of the atoms in the feedback region, around the 
estimated atom number N c : 

P 2 /(2mN w ) = P*/(2mN e )(l-5N w /N e +. . .), (6) 

where SN W = N w - N c . Following Ref. Q] we find that 
only the zeroth- and first-order terms of this expansion appear 
in the energy change AE: the expectation values of higher- 
order corrections exactly cancel each other. The leading term 
represents the cooling of the system by feedback. For this 
term, the energy removed is — (P 2 ) - / (2mN c ) . With m be- 
ing the atomic mass, N c m is the estimated total mass of the 
atoms in the feedback region, and P w is their total momen- 
tum. Thus this is the negative (estimated) kinetic energy of 
the center-of-mass of the atoms in the feedback region. Ac- 
cording to 0, this term contains atom-atom correlations of 
the form (0'(ri)$(r2)0'(r3)$(r4)) — a clear indication that 
in the quantum regime stochastic cooling cannot be described 
as a single-atom problem. 

The first-order correction in Eq. (jfji is found to give rise to a 
heating contribution to AE of the form (5N w P 2 )-/(2mN 2 ). 
From Eq. (jfji it can be seen, that this heating is due to non- 
optimal shifts of total momentum produced by atom-number 
fluctuations around the estimated value N c . Another way to 
see this is to consider a many-atom quantum state after a 
perfect measurement (a = 0) of P w with outcome P. It shall 
also be an eigenstate of TV^, with TV atoms in the feedback re- 
gion, i.e., I 'J) = \P,N). After the momentum shift the state 
becomes U(P) \P, N) = \P(1-N/N e ), N). That is, the mo- 
mentum will be shifted to zero only if N c = N. Given that 
the system in general is in a state of imprecise atom number, 
choosing iV c = (N w ) will only produce the correct momentum 
shift on average. Therefore, atom-number fluctuations in the 
feedback region are transferred into momentum fluctuations, 
that lead to the observed residual kinetic energy. This poses 
a fundamental limit to the perfect operation of the feedback 
loop and thus to stochastic cooling 1 13]. 

Moreover, since a measurement is involved in the con- 
trol loop, AE contains also measurement-induced heating 
terms: Since the total momentum is measured with resolu- 
tion cr, this fluctuations of the total momentum leads to a 
residual kinetic energy of the center of mass of the atoms 
in the feedback region. Together with the correction fac- 
tor (N w )-/N c that accounts for the difference between es- 
timated and average atom number, the resulting heating term 
is a 2 1 (2mN e ) x {N w )_/N c . Clearly a measurement of to- 
tal momentum with resolution a induces back-action noise 
oc cr -1 into the center of mass of the measured atoms. This 
noise leads to an increase in potential energy (i.e. heating) of 
the center of mass in the harmonic trap, which is found to be 
mN e ujQ / (8a 2 ) x (N w )_/N e , Given these contributions, the 
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FIG. 2: Energy change in units of Eq versus temperature. iVtot = 
10 , size of feedback region in x, y directions is Ago, locations are: 
centered (solid), shifted out of center by 5 x Ago (dashed). 
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FIG. 3: Temperature T (in /iK) versus number of feedback opera- 
tions n. Parameters as in Fig. 2, ujq/(2ti) = 500Hz. Solid (dashed) 
curve corresponds to centered (non-centered) feedback region, doted 
line indicates To. 



total change of energy due to a single feedback operation is 
then 

AE= mN c u* \ (N w )- (6N W PZ)- 

\2mN c 8ct 2 J N c 2m A c 2 2mN c ' 

(7) 

To achieve optimal cooling AE should be as negative as 
possible which is achieved by minimizing the heating terms. 
The optimal measurement resolution, cr op t = Apo^/W^, mini- 
mizes the measurement-induced heating [first term in Eq. 0] 
to ^Eq(N w ) _ /Nn, where Eq is the ground-state energy of 
a single atom II 211 . Choosing this value, the squared mea- 
surement resolution per atom equals the squared ground-state 
momentum uncertainty of a single atom, Apo- Since the num- 
ber of atoms in the feedback region changes with temperature, 
the estimate N c should have a temperature dependence. Thus 
<7 op t should be adapted during the cooling process, if possi- 
ble, to provide maximum cooling. For further optimisation the 
size and location of the feedback region in the xy-plane (cf. 
Fig. is of major importance. Here we consider only partic- 
ular cases; a more detailed study of the optimization will be 
presented elsewhere. In particular, we focus on temperatures 
near the condensation point, to demonstrate that in principle 
Bose-Einstein condensation can be reached with stochastic 
cooling. 

We have calculated the expectation values in Eq. Q using 
the grand-canonical ensemble for a system at temperature T 
and with average total number of atoms N tot . In this way 
we obtain the dependence of AE on the temperature for fixed 
iV to t, as shown in Fig. [2] for N tot = 10 6 atoms. If the feed- 
back region is centered with respect to the trap potential, it 
has substantial overlap with the condensate wavefunction. In 
this case (solid curve), when the temperature is gradually de- 
creased, there is a sudden change from cooling to heating at 
the condensation temperature To. This effect is due to dra- 
matic atom-number fluctuations of the condensate fraction at 
temperatures close to the phase transition 11411 . These fluctu- 



ations lead to a substantial heating due to the second term in 
Eq. 0. 

This problem can be avoided by choosing a feedback re- 
gion that has no spatial overlap with the nascent condensate 
and is thus not affected by its large atom-number fluctuations. 
In this case the energy removed per step gradually diminishes 
below T , though cooling still takes place (dashed curve). On 
the other hand, the energy subtraction per step at T ^ T is 
slightly smaller than for a centered feedback region. An ad- 
vantageous strategy is therefore to gradually move the feed- 
back region out of the trap center when approaching the con- 
densation temperature. 

Our present results lead us to the conclusion that neither 
quantum measurement effects nor the lack of knowledge of 
the precise number of atoms in the feedback region in princi- 
ple prevent Bose-Einstein condensation by stochastic cooling. 
However, we have said nothing about the speed at which con- 
densation may be reached. A conclusive answer on this issue 
requires that the cooling rate be dynamically calculated using 
the method described in Ref . 1 1 1 ] . For now we proceed with 
an equilibrium approach that represents a worst-case scenario: 
Starting at temperature Ti we calculate the energy subtraction 
AE(Ti), and thus a new average energy. From this energy a 
new temperature T i+1 is obtained by assuming the equilibra- 
tion of the system during the free evolution after the feedback. 
That is, we assume that due to collisions the atoms exchange 
energy and re-establish a new equilibrium state. Iterating this 
calculation many times we obtain the dependence of tem- 
perature on the number of feedback operations, as depicted 
in Fig. [3] Starting with a temperature of 7.6 /iK, condensa- 
tion is reached after about 10 7 feedback operations. Again a 
feedback region outside the trap center is advantageous when 
reaching T (dashed curve). 

In a realistic experiment stochastic cooling is performed 
without waiting for re-equilibration between feedback oper- 
ations 1 15]. Then non-vanishing coherent amplitudes of the 
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momentum (P w ) will appear depending on the oscillation 
phase, i.e., the (randomly chosen) time of the subsequent feed- 
back operation. Since (P 2 ) = ((AP W ) 2 )+(P W ) 2 , the feedback 
then not only reduces momentum fluctuations ((AP^,) 2 ), as 
for a thermal state, but also subtracts the energy due to the 
non-equilibrium coherent amplitudes (P w ) 2 [cf. Eq. Q]. The 
latter leads to additional energy subtractions as compared to 
our equilibrium approach. One can thus expect a much faster 
cooling process, so that our results (cf. Fig. |5J represent the 
upper limit of the required number of feedback operations. 

Our calculations have shown that increasing the size of the 
feedback region further increases the cooling efficiency. Un- 
fortunately, for enlarged feedback regions an increased num- 
ber of trap levels is required, which presently runs into lim- 
itations of our numerics. For finding optimal strategies and 
parameter ranges of stochastic cooling we intend to imple- 
ment in the near future a dynamical calculation 1 11]. Then 
the dependence of the cooling rate on size and location of the 
feedback region can be studied in full detail. 

It is worth noting, that the considered feedback regions 
could be easily realized in experiment by application of op- 
tical fields [ 10]. The ground-state position variance of sodium 
atoms is approximately Ago = 0.6/im for a trap of frequency 
v = 500Hz. Beam waists of externally applied laser fields, that 
implement measurement and shift, can be chosen in a wide 
range limited only by optical wavelengths, which are much 
smaller than Ago- Thus sizes of the feedback region much 
smaller and much larger than Ago could be realized. 

Note also, that the geometry used here is different from 
that in Ref. 1 10]. Raizen et al. considered a one-dimensional 
model where the feedback region was restricted in that same 
direction. They concluded that finer spatial resolution of the 
feedback region leads to increasing cooling efficiencies. In 
our 3D geometry that scenario would correspond to a restric- 
tion of the feedback region also in z direction, say to a size 
Az. Then the measurement can resolve momenta only within 
a resolution that is enlarged by l/(2Az). This however, may 
decrease the amount of subtracted energy. The latter effect is 
absent in the classical calculation as performed in Ref. iTcfc . 
For now we cannot confirm the predicted increase of cool- 
ing efficiency with increased spatial resolution in z direction, 
since such a scenario would lead to non-trivial modifications 
in Eq. 0. 

In conclusion by using a quantum-field theoretical ap- 
proach we have derived the energy change of feedback opera- 
tions of stochastic cooling of trapped bosonic atoms. Besides 
the heating due to the quantum measurement and the sought 
subtraction of kinetic energy, we have shown that stochas- 
tic cooling is strongly governed by a noise term that is due 
to atom-number fluctuations in the feedback region. This ef- 
fect becomes dominant at the condensation point where atom- 
number fluctuations are large. This detrimental heating effect 
can be ameliorated by the choice of the feedback region. It 
has been further shown that condensation temperatures can in 
principle be reached. Our results are based on an equilibrium 
approach and higher cooling efficiencies are expected for a 



non-equilibrium dynamics as realizable in experiment. Fu- 
ture investigations with fully dynamically solutions will pro- 
vide further insight into stochastic cooling with respect to op- 
timization of cooling rates. 
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